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After taking into account both the pion and sigma meson exchange potential, we have performed 
a dynamical calculation of the D'^D*'^ system. The a meson exchange potential is repulsive from 
heavy quark symmetry and numerically important for a loosely bound system. Our analysis disfavors 
the interpretation of X(3872) as a loosely bound molecular state if we use the experimental D'Dtt 
coupUng constant g = 0.59 and a reasonable cutoff around 1 GeV, which is the typical hadronic 
scale. Bound state solutions with negative eigenvalues for the DD* system exist only with either 
a very large coupling constant (two times of the experimental value) or a large cutoff (A ~ 6 GeV 
or /3 ~ 6 GeV^). In contrast, there probably exists a loosely bound S-wave BB* molecular state. 
Once produced, such a molecular state would be rather stable since its dominant decay mode is the 
radiative decay through B* — > By. Experimental search of these states will be very interesting. 

PACS numbers: 12.39.Pn, 12.40.Yx, 13. 75. Lb 



I. INTRODUCTION 

Since the observation of the charmonium-like state 
X{3872) in the J/i/itt+tt" channel by Belle collaboration 
in 2003 [H, X(3872) has been confirmed by CDF 0, DO 
Q and Babar collaborations In the past three years, 
there have accumulated abundant experimental informa- 
tion of X(3872), which is collected in TablcU 



Mass 
(MeV) 



Width 
— 



Decay channels 



Branching fractions 



X(3872) 



3872.0 ±0.6 ±0.5 [IJ 
3871.3 ±0.7 ±0.4 [2] 
3871.8 ±3.1 ±3.0 [3] 
3873.4 ± 1.4 [4] 

3875.4 ± 0.7^2.? [5] 
3875.6 ± 0.7li1 [6] 



< 2.3 MeV [1] 



1++/2-+ [7, 8] 



X(3872) ^ J/tpTv+Tv- [1,2,^, J 
X(3872) ^ yJ/ip,u}J/tp'^'M 
A" (3872) ^ pJ/jP [llj 
A(3872) ^ D°D°n° [5] 
A(3872) ^ D°D*° + h.c. [6] 



-B-R|X(3872)-.7J/;/)] _ 
SH[X(3872)-.7r+ir- J/i>] ^ 
BR[X(3872)^-fJ/^] 

Sfl[X(3872)-.7r+7r- J/V] 
BR[X(3872)^D° D°tt"] 



BR\X(3872)- 



0.14 ±0.05 [9] 
= 0.25 [10] 
= 9.4iij [5] 



TABLE I: A review of the experimental status of A(3872). 

Quark model calculation indicates that a 1?P\ cc state 
X'ci lies 50 - 200 MeV above X(3872). Moreover a 
charmonium state with isospin / = docs not decay 
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into J/ipP easily. Thus there is some difficulty of the 
charmonium assignment of X(3872). The possible theo- 
retical explanations of X(3872) include a molecule state 
[13, [13, Q Ed, a 1_++ cusp [l3, the S-wave threshold 
effect due to the D°D°* threshold [H], a hybrid char- 
monium [l3| , a diquark anti-diquark bound state [l^] , a 
tctraquark state [2l| and a dynamically generated reso- 
nance [2^ . 

Among these theoretical schemes, the molecule picture 
is the most popular one due to the following reasons. The 
molecular picture naturally explains both the proximity 
of a: (3872) to the D°D*° threshold and the isospin vi- 
olating J/ipp decay mode. It predicted the decay width 
of the J / ipTT^ TT~ n'^ mode to be comparable with that of 
J/'ij^p, which was confirmed by Belle collaboration [3]. 
Within the same picture, Braaten and Kusunoki pre- 
dicted that the branching ratio of X(3872)A'° is 
suppressed by more than one order of magnitude com- 
pared to that of B+ -> X{3872)K+ 

Later both Belle and Babar collaborations observed 
the radiative decay mode. Belle's measurement found [3| 



BR[X(3872) -fJ/tP] 
BR[X{3872) J/V'TT+TT-] 

while Babar collaboration got [l3| 

BR[X{3872) yJ/^j] 



BR[X{3872) ^ J/^m+ir- 



0.14 ±0.05 



0.25 , 



(1) 



(2) 



which are against the prediction from the molecular pic- 
ture 7 X 10~3. 

Recently Belle collaboration measured the ratio [3] 



BR[X{3872) D^D^TT^] 
BR[X{3872) TT+ir-J/^j] 



— q 4+3.6 



(3) 



which is much larger than the theoretical value 0.054 
from the molecular assumption. From Ref. [3], one can 
also extract 



BR[B" X{3872)K"] 
BR[B+ X{3872)K+] 



1.62 



(4) 



2 



which is also much larger than the molecule prediction. 

Up to now, several groups carried out the dynamical 
study of the molecular assignment of X(3872). Swanson 
proposed that X{3872) was mainly a D^D*^ molecule 
bound by both the pion exchange and quark exchange 
(l5l |. To obtain the potential between D°D*^ through 
exchanging single pion, he followed the method proposed 
by Tornqvist |24i |. The formalism is based on a micro- 
scopic quark-pion interaction. Swanson indicated that 
one pion exchange alone can not bind D and D* . He 
also included the short-range quark-gluon force • 

In Ref. [ill, Wong studied the DD* system in the 
quark model in terms of a four-body non-relativistic 
Hamiltonian with pairwise effective interactions. This 
framework is similar to the consideration of adding short- 
range quark-gluon force in Swanson's paper [l5|. The 
author found an S-wave DD* molecule with the binding 
energy ~ 7.53 MeV. In Refs. [H, [H, [13, [H, [H, , fur- 
ther investigations basing on the molecular assumption 
are carried out. 

With the obtained one pion exchange potential 
(OPEP) by using the effective Lagrangian, Suzuki argued 
that X(3872) is not a molecular state of D°D*° + D°D*° 
[3^ . which contradicts Swanson and Wong's conclusion. 
Instead, X(3872) may have a dominant cc component 
with some admixture of D°D*^ + D°D*^ [3l|, [H, IH . 

In order to further clarify the underlying structure of 
X(3872), we shall carry out a systematic dynamical study 
of the molecular picture in this work. It's important to 
note that the one pion exchange potential alone does not 
bind the proton and neutron pair into the deuteron in 
nuclear physics. In fact, the strong attractive force in 
the intermediate range has to be introduced in order to 
bind the deuteron, which is modelled by the sigma meson 
exchange potential elegantly. We shall explore whether 
the similar mechanism plays an important role in the case 
of X(3872). 

This work is organized as follows. After the intro- 
duction, we give a concise review of the molecular pic- 
ture. In Section IlIII we present the flavor wave function 
of X(3772), effective Lagrangian and coupling constants 
relevant to the derivation of the tt and a exchange po- 
tentials. In Section IIV| we illustrate the procedure to 
obtain the potentials and give their expressions. Then 
we present the numerical results in Section |V| and IVII 
The last section is the summary and discussion. 



II. REVIEW OF MOLECULAR PICTURE 

In the study of hadron spectroscopy, some states are 
difficult to be accommodated in the conventional qq and 
qqq framework. These states are considered good candi- 
dates of hadrons beyond the conventional valence quark 
model. The possible assignments include the glueball, 
hybrid state and mutiquark state etc. Among them, the 
molecular state is very attractive. 

In the past thirty years, theorists have been study- 



ing whether two charmed mesons can be bound into the 
molecular state because the presence of the heavy quarks 
lowers the kinetic energy while the interaction between 
two light quarks could still provide strong enough attrac- 
tion. Voloshin and Okun studied the interaction between 
a pair of charmed mesons and proposed the possibilities 
of the molecular states involving charmed quarks [s^ . 
dc Rujula, Gcorgi and Glashow once suggested -0(4040) 
as a D*D* molecular state [35|. Tornqvist studied the 
possible deuteron-like two-meson bound states such as 
DD* and D*D* using the quark-pion interaction model 
[23 |. Dubynskiy and Voloshin indicated that there exists 
a possible new resonance at the D*D* threshold [H, [IJ- 
Besides the above systems, Weinstein and Isgur studied 
whether the scalar resonances /o(980) and ao(980) are 
molecular states composed of a pair of KK mesons [s^ . 

In the past several years, the experimental observations 
of so many X, Y and Z states stimulated the study of ex- 
otic states greatly. For example, AT (3872) is proposed to 
be a good candidate of the DD* molecule state by many 

groups [12, [13, [3, [H, [3, which is also the topic of the 

present work. Liu, Zeng and Li suggested Y(4260) as the 
Xcp'^ molecule assignment and predicted its possible de- 
cays modes [11] . Yuan, Wang and Mo proposed F(4260) 
being a Xa^ molecule [13]. The baryonium possibility 
was also suggested by Qiao [4ll |. 

Recently Belle collaboration observed a charged state 
Z+(4430) in t/j'tt channel [i^. This new enhancement 
immediately triggered the molecular speculation. In fact 
several g rou p s suggested Z+(4430) as a DiD* molecular 
state [43, [ill. In our previous work [4^, we have carried 
out the first dynamical study of ^+(4430). Later, we 
performed a detailed study of this state in the molecular 
picture [i^. A short review of the current theoretical 
status of Z+ (4430) [H, [H, [H, [13, [HI, [13, [H, [13, 

[55! was also given in Ref. j45j |. 



III. FLAVOR WAVE FUNCTION, EFFECTIVE 
LAGRANGIAN AND COUPLING CONSTANTS 

In the following, we will study whether X(3872) is a 
bound state of the DD* meson pair. Before deriving 
the meson exchange potential, we first briefiy discuss the 
convention of the flavor wave function of the molecular 
state X(3872). In the previous hterature [13, [H, [13, [H, 
[l6j . it was defined as 



|X(3872)) 



1 



7=)*0\ 



WD 



c\D*"D°) 



(5) 



with c = +1. However, this definition does not refiect 
the positive C-parity of X(3872) naturally [111- Accord- 
ing to the same approach in our previous paper [45| . we 
reanalyze the fiavor wave function of X{3872). 

The interpolating current of X(3872) corresponding to 
Eq. ([5]) in the quantum field theory reads 



J 



I 



Jf(3872) 



V2 



( Ji + CJ2) 



(6) 



3 



with 



Ji 



(u"75c")(cS^w''), J2 = (c''75M°)(mS^c''), 



where a, 6 denotes the color indices. Under the charge 
conjugate transformation, one gets 

CJiC-^ = -J2 and CJ2C-^ = -Ji. 

We want to emphasize that there exists no arbitrary 
phase because the charm and anti-charm quark and the 
up and anti-up quark appear simultaneously. Therefore 
we obtain 



estimate the values of the coupling constant g^, we com- 
pare the Lagrangian with that in Ref. [s^ and get 



9<y = 



2V6 



(12) 



with = 3.73. Unlike the case of Z+(4430) [4a|, it 
is unnecessary to care about the phases of the coupling 
constants in the present case. We will turn to this point 
later. 



CJ 



X(3872 



1 



-.{-J2-CJI). 



IV. THE DERIVATION OF THE ONE PION 
AND SIGMA EXCHANGE POTENTIAL 



Because the charge parity of X(3872) is -1-1, we have 
c = — 1. In other words, the natural definition of the 
flavor wave function of X(3872) should be 



1 



|X(3872)) = 



\D°D*") - \D*°D") 



(7) 



In this work, we mainly discuss whether the S-wave D'^ 
{D°) and D*° {D*°) molecular state can be formed by 
exchanging the tt and a meson. We need the effective chi- 
ral Lagrangian in the chiral and heavy quark dual limits 



with 



igTT[Hb4ba75Ha] + gaTvlHaH] 
1+ 



Ha = 



and the axial vector field is defined as 
with ^ = exp{iM/.U), .U = 132 MeV and 



(8) 



(9) 



To derive the effective potential, we follow the same 
procedure in Ref. [i^ . Firstly we derive the elastic scat- 
tering amplitudes of both the direct process and crossed 
channel. Secondly, we get the potential in the momen- 
tum space for a special component (e.g. Jz = 0) with 
the Breit approximation. Then we average the poten- 
tial in the momentum space. Finally we make Fourier 
transformation to derive the potential in the coordinate 
space. 

In the present case, the parity and angular momentum 
conservation ensures that the tt exchange occurs only in 
the crossed channel while the a exchange only in the 
direct channel (see Fig. [T]). The zeroth component of 
exchange meson momentum is qo ~ Mi — Mf. For the 
direct scattering diagram, Mij denotes the mass of 
Thus we can approximately take qo = and q- = — q". 

However, qo could not be ignored because Mi and Mf 
denote respectively the masses of D° and D*° for the 
crossed diagram, qo ~ M^-o — Aljja is larger than pion 
mass TOtt which indicates that the exchanged pion can be 
on-shell. In this case, one can deal with the potential in 
the coordinate space by the principal integration as in 
Eq. ([H]) below. 



M = 



I 



+ 



_rL 

\/2 V6 

K- 



7r+ if+ 



(10) 



In Ref. [56| , the coupling constant g = 0.75 was esti- 
mated roughly within the quark model. A different set 
of coupling constants can be found in Ref. [sH]. With 
our notation, g ~ 0.6 In fact, the coupling constant 
g was studied using many theoretical approaches such as 
QCD sum rules jsol. [60. 61, 62]. Despite so many theo- 
retical estimates of the coupling constant g, we use the 
value 



<7 = 0.59 ±0.07 ±0.01 



(11) 



in this work. The above value was extracted by fitting 
the precise experimental width of D* [63| . In order to 



r," . 



FIG. 1: The scattering of D" — D*" by exchanging the tt and 
CT mesons. 



We use the following definitions in the potentials after 
Fourier transformation 



1 C^qr 



q2 + ml 



(27r)3 



q2 _ ^2 



(2^ 



(13) 
(14) 



4 



Writing them explicitly, we have 



1 



y^(r) = -(5(r) - — cos(/^r), 



(15) 



where fj. = i/gg — m^. Except the relative sign, y,r(r) 
is similar to the expression derived in Ref. [sij by using 
the polarization vectors e^^ = -^(0,±l,i,0) and e° ~ 

(0,0,0,-1) M- 

With the convention of the X(3872) flavor wave func- 
tion in Eq. ([7]) , the potential in the study of the molecular 
picture finally reads as 



y(r)=5^y.(r) + |-F.(r). 



(16) 



Here the sign between one sigma exchange potential 
(OSEP) and OPEP is determined by the relative sign 
of \D°D*°) and \D*°D°) in the wave hmction in Eq. 1^. 

It's important to note that the signs in the potential 
are completely fixed. The heavy quark spin-flavor sym- 
metry ensures that the D and D* mesons possess the 
same coupling constants. The resulting potential in Eq. 
(fTB)) does not change with the phases of coupling con- 
stants. 

Especially, wc find that a exchange potential is repul- 
sive, which differs from that in the nuclear forces. Be- 
cause of this unique feature, one just needs to study 
whether the one-pion exchange can bind D and D* 
mesons to form X(3872). Only when the answer is posi- 
tive, should we consider the effect from the a exchange. 

We note that the potential in Eq. p6|) is derived with 
the implicit assumption that all the mesons are point-like 
particles. Such an assumption is not fully reasonable due 
to the structure effect in every interaction vertex depicted 
in Fig. [TJ Thus in the following we will introduce the 
cutoff to regulate the potential and further study whether 
it is possible to find a loosely bound molecular state using 
the realistic potential. 

We will modify the potential through two approaches: 
(1) considering the form factor (FF) contribution; (2) 
smearing the potential. Although these two approaches 
look different, they are essentially the same, i.e. imposing 
a short-distance cutoff to improve the singularity of the 
effective potential. 



A. Introducing form factors in the potential 

Before making a Fourier transformation, we introduce 
a form factor in the interaction vertex to compensate the 
off-shell effects of the exchanged mesons. The adopted 
FF is of the monopole type [2J, [g^] 



A^ — m' 
A2-g2 



(17) 




FIG. 2: The regulated potentials related with X{3872) in the 
case of FF. The solid line corresponds to OPEP. The long- 
dash line comes from one a exchange, and the short-dash 
line is the total effective potential, g = 0.59, Qa = 0.76 and 
A = 1.0 GeV are used. 



where A ^ 1 GeV denotes a phenomenological cutoff, m 
and q are the mass and the four-momentum of the ex- 
changed meson respectively. As 0, FF becomes 
a constant. With A ^ m, it approaches unity. In 
other words, as the distance is infinitely large, the vertex 
looks like a perfect point. So the form factor is simply 
unity. On the other hand, as q^ oo, the form fac- 
tor approaches to zero. In this situation, as the distance 
becomes very small, the inner structure (quark, gluon 
degrees of freedom) would manifest itself and the whole 
picture of hadron interaction is no longer valid. 

The explicit expressions of the modified potentials are 



1 

47rr 



-Ar 



-Ar 



8nA 



(18) 



Y^r) = -|i[cos(H-e-"1-^e-"^ (19) 



where 77 = yA^ 



~ — rri^ and 



v/A2 - ql. Note we use the same A for tt and a exchange. 
As an example, we have plotted the above regulated po- 
tential in Fig. [51 



B. Regulating the potential with the smearing 
technique 

The potential can be written as 

V{t) = j V{y')6{v - r')dv' . (20) 
To smear the potential, we employ the replacement 



5{y-v') 



3/2 



(21) 



5 



— j[ 

— - o 

j[+o 




r[GeV"'] 



To solve the Schrodinger equation, one needs the fol- 
lowing parameters: m^r = 134.98 MeV, nig- = 600 MeV, 

= 132 MeV, mo- = 2006.7 MeV, m^o 1864.6 MeV 
[67|. In this seetion, we first eonsider whether the one 
pion exchange interaction alone ean bind DD*. 

Now we explore at what condition D and D* can form 
a bound state through one pion exchange interaction with 
two approaches. Our procedure to collect the numerical 
values is: (1) we fix the coupling constant g = 0.59 and 
vary the cutoff (A or /3) from a small value until we find 
a solution with a binding energy less than 5 MeV; and 
(2) we increase g to several larger numbers and tune the 
cutoff until a solution with a binding energy less than 5 
MeV is found. 



A. Results for the case of FF 



FIG. 3: The regulated potentials related with X{3872) in the 
case of smearing. The solid line corresponds to OPEP. The 
long-dash line comes from one a exchange, and the short-dash 
line is the total effective potential, g = 0.59, = 0.76 and 
/3 = 1.0 GeV^ are used. 

which was suggested by Isgur in Ref. [11]. As (3 goes 
to infinity, the right-hand-side of the above expression 
becomes the delta function. Typical values of ^/P are 
~ 1 GeV, corresponding to the short range cutoff. 
I.e., the short-distance structure is indiscriminate. On 
the other hand, /3 should not be very small to describe a 
system with internal structure. 
We obtain the smearing potential 

V{r) 

smearing 

gl _pr'\ ''"%,f'->' . f m^- Wr \ 

= e ' e V erfc — — 

87rr L V 2V^ / 

V 2V^ )\ 6/2 v^y 

- ^ , e 15 erf ' ^ + c.c. . 22) 

Here erf(x) and erfc(x) denote the error function and 
complementary error function respectively while c.c. de- 
notes complex conjugate. An illustrative example of the 
smeared potential is presented in Fig. [3] 



V. NUMERICAL RESULTS FROM THE ONE 
PION EXCHANGE INTERACTION ALONE 

In order to find whether there is a bound state in DD* 
system, we solve the radial Schrodinger equation with the 
help of MATSLISE [e^, which is a graphical MATLAB 
package for the numerical solution of Sturm-Liouville and 
Schrodinger equations. A bound system has at least one 
negative eigenvalue. 



If the coupling constant g is fixed to be the experi- 
mental value g = 0.59, the possible bound state solu- 
tion with a negative eigenvalue can only be found when 
A > 5.6 GeV. The larger the cutoff A is, the closer the 
regulated potential is to the delta function, hence the 
larger the binding energy. The binding energy is very 
sensitive to A. This result is consistent with the behav- 
ior that F{q^) 1 when A oo. It's known that the 
three-dimensional — (5(r) function alone does not gener- 
ate a bound state. The requirement A > 5.6 GeV is 
much much larger than the commonly used reasonable 
value ~ 1.0 GeV. In other words, the one pion exchange 
potential alone does NOT bind the D'^D*'^ pair into a 
molecular state with the physical values of g and A! This 
is our first important observation. 





A (GcV) Eo (MeV) 


rrms (fm) 


rmax (fm) 


g = 0.59 


5.7 


-0.3 


5.8 


0.2 


5.8 


-2.1 


2.2 


0.2 


3 = 0.7 


4.1 


-0.8 


3.7 


0.2 


4.2 


-3.2 


1.8 


0.2 




3.1 


-0.1 


8.7 


0.4 


3 = 0.8 


3.2 


-1.6 


2.6 


0.3 




3.3 


-4.9 


1.5 


0.2 


3 = 0.9 


2.5 


-0.6 


4.2 


0.4 


2.6 


-2.9 


2.0 


0.3 


3 = 1.0 


2.0 


-0.2 


7.2 


0.5 


2.1 


-1.8 


2.5 


0.4 



TABLE 11: Solutions for various g and A in the case of FF 
with OPEP. Lowest eigenvalues between -5.0 MeV and -0.1 
MeV are selected. 

We consider only the solutions with the eigenvalues be- 
tween -0.1 MeV and -5.0 MeV corresponding A = 5.7 and 
A = 5.8. To understand the solutions more clearly, we 
present the numerical results in Table |TT1 Eo is the low- 
est eigenvalue of the system, rmis is the root-mean-square 
radius, and rmax is the radius corresponding to the max- 
imum of the wave function In Fig. 3] and [51 we 
present the radial wave functions R{r) and x(r) = rR{r) 
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A.5.7 GeV 

- - A.5.8 GeV 




A=5.7 GeV 

- - A=5.8 GeV 



FIG. 4: The radial wave functions R(r) corresponding to A : 
5.7 GeV and A = 5.8 GeV with g = 0.59. 



respectively. According to the figures, as A increases, the 
probability for a bound state appearing near the origin 
becomes larger. The large value of rmis indicates this 
possible bound state is very extended, which can be il- 
lustrated with the figures. 

Secondly, we enlarge g arbitrarily until g = 1.0 and 
perform a similar evaluation. The results are also pre- 
sented in Table |lTl When g becomes larger, the crit- 
ical point for A to generate a DD* bound state be- 
comes smaller. With a reasonable cutoff A ~ 1.0 GeV, a 
bound state exists only when the coupling is very strong 
{g > 1.0), which is nearly two times of the experimental 
value. The wave functions corresponding to the solutions 
in Tabic [TTl have similar shapes with those in Figs. [4] and 

El 

Now we come back to discuss the partner state of 
X(3872). We denote it as X. The C parity of X is 
negative. 



\X) = 



1 

V2 



f=i*0\ 



\D*"D°) 



(23) 



With this convention, the signs in the OPEP are reversed 
while the sigma meson exchange is still repulsive. There- 
fore the attractive force is much weaker. We find that the 
potential is not attractive enough to bind D and D* even 
with g = 1.0. If we arbitrarily use g = 5.0 and A = 1.0 
GeV, one finds a negative eigenvalue about -0.1 McV. 
The value is not sensitive to A. In this case, rmis ~ 19 



fm and 



is about 14 fm. From these values, one 



concludes that this convention does not lead to a DD* 
bound state with the realistic coupling constant. It is not 
difficult to understand the results with the potential in 
Eq. 1151 The part which could provide some attraction 
is ^^2^^. Since ^ = 0.044 GcV is smah, a possi- 
ble bound state exists only if g is a very large number. 
The consideration of FF improves mainly the behavior 




r[GeV-'] 

FIG. 5: The function x{^) ~ rR{r) corresponding to A = 5.7 
GeV and A — 5.8 GeV with g = 0.59. The lower diagram 
shows the behavior in short range. 



of the most singular part. Thus the binding energy is 
insensitive to the cutoff. 

From the above analysis, we conclude that DD* in- 
teraction through one pion exchange is not attractive 
enough to form a bound state with g = 0.59 and A 1.0 
GeV. 



B. Results for the case of smearing 

In the case of the smeared potential, one fails to 
find a bound state solution with negative eigenvalue for 
(3 < 5.3 GeV^ if we fix 5 = 0.59. The binding energy is 
very sensitive to and increases with f3. With a reasonable 
cutoff /? ~ 1 GeV^ , there exists no loosely bound molec- 
ular state using the realistic coupling constant g ~ 0.59. 

When we vary g from 0.59 to 1.0 and select the solu- 
tions with -5.0 MeV < £^0 < -0.1 MeV, we obtain the 
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P (GeV") 


Eq (MeV) 


rrms (fm) 


r-max (fm) 






5.5 


-0.3 


5.8 


0.2 




= 0.59 


5.6 


-1.0 


3.3 


0.2 


g 


5.7 


-2.0 


2.3 


0.1 






5.8 


-3.4 


1.8 


0.1 






2.8 


-0.3 


5.5 


0.2 


g 


= 0.7 


2.9 


-1.4 


2.7 


0.2 






3.0 


-3.3 


1.8 


0.2 




= 0.8 


1.7 


-0.9 


3.4 


0.3 


g 


1.8 


-3.1 


1.9 


0.2 




= 0.9 


1.1 


-1.4 


2.8 


0.3 


g 


1.2 


-4.9 


1.5 


0.3 




= 1.0 


0.7 


-0.5 


4.4 


0.4 


g 


0.8 


-3.9 


1.7 


0.3 



TABLE III; Solutions for various g and (5 in the case of smear- 
ing with OPEP. Lowest eigenvalues between -5.0 MeV and 
-0.1 MeV are selected. 

results in Tabic IIIII One gets similar conclusion as in 
the form factor case. The critical point for (3 to generate 
a bound state is lowered as g becomes larger. For ex- 
ample, a bound state can be obtained with g ~ 0.9 and 
/3 ~ 1.0 GeV^. The shapes of the wave functions corre- 
sponding to these solutions are also similar to those in 
Figs. HandEl 

As in the form factor case, if the flavor wave function 
is used, no bound states can be found with g = 1.0. 
If 5 = 5.0, a bound state exists and the eigenvalue is 
insensitive to the cutoff. The numerical results are very 
close to those in the form factor case, which also indicates 
the insensitivity of the results to the cutoff. Therefore, 
it is also difficult to find a DD* bound state by one pion 
exchange interaction with the realistic coupling constant 
g — 0.59 in the smearing case. 

From the above analysis within two approaches, 
we find that the molecular interpretation of A'(3872) 
through one pion exchange interaction may be problem- 
atic. The regulated OPEP may generate bound states ei- 
ther with an unphysically large coupling constant i? > 1.0 
or an un-reasonably large cutoff. The bound state solu- 
tion with the realistic coupling constant does not exist 
if the value of the cutoff is around 1 GeV. The two ap- 
proaches agrees with each other and lead to the same con- 
clusion. As a by-product, we point out that our sign con- 
vention for the flavor wave function of X(3872) is much 
more helpful to form a bound state than the old conven- 
tion used in the literature. 



VI. NUMERICAL RESULTS WITH BOTH THE 
PION AND SIGMA MESON EXCHANGE 
INTERACTION 

Now we move on to include the one a exchange inter- 
action. The (7 contribution reinforces the above conclu- 
sion in the previous section due to the repulsive nature 
of OSEP. We will study carefully the variation of the nu- 
merical results and see how much it affects the conclusion 



when OSEP is considered. The procedure is similar to 
the OPEP case. 



A. Results for the case of FF 

We first take a look at the potentials plotted in Fig. [51 
The curves are obtained with g = 0.59, g^ = 0.76, and 
A = 1.0 GeV. From this figure, one notes that OSEP is 
small compared with OPEP. Thus one expects one sigma 
exchange interaction has small contributions to the bind- 
ing energy. However, since a very loosely molecular state 
is expected, a small variation of the potential may lead 
to relatively big change of the eigenvalue. 

By adding OSEP in the Schrodinger equation, one gets 
numerical solutions listed in Tabic [TVl We only use the 
coupling constant g^ = 0.76 to illustrate the results. 
Again, we chose the solutions with —5.0 MeV < Eq < 
-0.1 McV. 





A (GeV) Eo (MeV) 


rrms (fm) 


rmax (fm) 


g = 0.59 


6.0 


-1.3 


2.8 


0.1 


6.1 


-4.9 


1.5 


0.1 


3 = 0.7 


4.3 


-1.1 


3.1 


0.2 


4.4 


-4.5 


1.5 


0.2 


3 = 0.8 


3.3 


-0.7 


3.8 


0.3 


3.4 


-3.7 


1.7 


0.2 


3 = 0.9 


2.6 


-0.4 


5.0 


0.3 


2.7 


-2.8 


2.0 


0.3 


3 = 1.0 


2.1 


-0.3 


5.9 


0.4 


2.2 


-2.4 


2.2 


0.3 



TABLE IV: Solutions for various g and A in the case of FF 
with total potential. Lowest eigenvalues between -5.0 MeV 
and -0.1 MeV are selected. Here g^ = 0.76 is used. 

By comparing the data in Tables |TT] and llVi one finds 
that many bound state solutions with negative eigenval- 
ues for certain pairs of g and A disappear after we include 
the repulsive sigma meson exchange force. Only three so- 
lutions survive with -5.0 MeV < £^0 < -0.1 MeV. But 
their binding energy decreases by at least 83%, which 
clearly indicates that the sigma exchange force are nu- 
merically very important for a loosely bound molecular 
state. 



B. Results for the case of smearing 

The smeared potentials is plotted in Fig[3l where we 
use g = 0.59, P = IGeV^ and g^ = 0.76. By using 
gcr = 0.76 and selecting solutions for Eq between -5.0 
MeV and -0.1 MeV, we get the results given in Table 
IVl Comparing data in this table with those in Table 
mil only two solutions (when g = 0.9, /3 = 1.2 GeV^ and 
g = 1.0, P = 0.8 GeV^) still satisfy our requirement. The 
binding energy decreases by at least 74%. 
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P (GeV") Eo 


(MeV) 


rrms (fm) 


r-max (fm) 






6.0 


-0.1 


8.6 


0.1 




= 0.59 


6.1 


-0.8 


3.7 


0.1 


g 


6.2 


-1.9 


2.4 


0.1 






6.3 


-3.5 


1.7 


0.1 






3.1 


-0.3 


5.5 


0.2 


9 


= 0.7 


3.2 


-1.6 


2.6 


0.2 






3.3 


-3.8 


1.7 


0.2 




= 0.8 


1.9 


-1.2 


3.0 


0.2 


9 


2.0 


-3.9 


1.7 


0.2 




= 0.9 


1.2 


-1.0 


3.2 


0.3 


9 


1.3 


-4.4 


1.6 


0.3 


9 


= 1.0 


0.8 


-1.0 


3.3 


0.3 



TABLE V: Solutions for various g and (5 in the case of smear- 
ing with total potential. Lowest eigenvalues between -5.0 MeV 
and -0.1 MeV are selected. Here — 0.76 is used. 



VII. NUMERICAL RESULTS FOR BB* SYSTEM 



Finally we apply the formalism to BB* system. 



1 



\B+B*-) - \B*+B- 



(24) 



Because of the heavier masses of the B mesons, the kine- 
matic term has relative small contribution. The pos- 
sibility of forming a bound state is larger than that 
in the DD* system. OSEP remains the same. But 
the expression of the OPEP is different now because 
= mB* — TUB < rriTr. Therefore the potential can 
be strictly derived and does not have an imaginary part. 
Now we have 



-6{r) + fie- 



where fiB = \J'm\- {q%)^. 

If a form factor is introduced before the Fourier trans- 
formation, this function becomes 



where as = -^A^ — {q%Y ^^'^ V — ^/A^ — m^. 

If the smearing technique is applied, this function is 
regulated as 



YAr) 



fj\ 3/2 2 

87rr 

MB - 2/3r 



4/3 



erfc 



rfc 



/ijj + 2/3r 



A.2.3 GeV 

- - A.2.4 GeV 




r [GeV"'] 



FIG. 6: The radial wave functions R(r) corresponding to A = 
2.3 GeV and A = 2.4 GeV with g = 0.59 for the BB* system. 





A (GeV) Eo (MeV) 


rrms (fm) 


rmax (fm) 


g = 0.59 


2.3 
2.4 


-0.9 
-2.8 


2.1 
1.2 


0.3 
0.3 


3 = 0.7 


1.7 


-0.8 


2.3 


0.4 


1.8 


-2.7 


1.3 


0.3 




1.3 


-0.1 


5.1 


0.6 


3 = 0.8 


1.4 


-1.4 


1.7 


0.4 




1.5 


-4.2 


1.1 


0.4 


3 = 0.9 


1.1 


-0.4 


3.2 


0.6 


1.2 


-2.3 


1.4 


0.4 


3 = 1.0 


1.0 


-1.5 


1.7 


0.5 


1.1 


-5.0 


1.0 


0.4 



TABLE VI: Solutions for various g and A in the case of FF for 
(25) the BB* system with OPEP. The lowest eigenvalues between 
-5.0 MeV and -0.1 MeV are selected. 



are the same as in the DD* case. With the same proce- 
dure as before, we obtain solutions in various cases. Re- 
sults from the one pion exchange interaction for the case 
(26) of FF (smearing) are presented in Tabic IVTl (jVIip . After 
considering effects from the one sigma exchange interac- 
tion, the results corresponding to the case of FF (smear- 
ing) are collected in Table IVIIII (|IXp . For comparison, 
we also present the radial wave function R(r) and x('") i^i 
Figs, inland [71 From these tables, it's very interesting to 
note that there probably exists a loosely bound S-wave 
BB* molecular state. Once produced, such a molecu- 
lar state would be rather stable since its dominant decay 
mode is the radiative decay through B* Bj. 



(27) VIII. SUMMARY AND DISCUSSIONS 



When performing numerical evaluations, mB* = 5325 
MeV and niB ~ 5279 MeV [b^]. For the coupling con- 
stants, we use the values in the heavy quark limit which 



In this work we have studied whether X(3872) is an 
S-wave DD* molecule state bound by the one pion and 
one sigma exchange interactions. Wc choose to work at 
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A=2.3 GeV 

- - A=2.4 GeV 




FIG. 7: The function xi''') ~ rR{r) corresponding to A = 2.3 
GeV and A = 2.4 GeV with g = 0.59 for the BB* system. 
The lower diagram shows the behavior in short range. 



the hadronic level and employ the effective Lagrangian 
incorporating both the heavy quark symmetry and chiral 
symmetry. We find the a meson exchange potential is 
repulsive and numerically important for a loosely bound 
system. 

Considering the internal structure and finite size of the 
hadrons, we have regulated the singular S function in the 
potential using both the form factor and smearing tech- 
nique. After solving the radial Schrodinger equation with 
regulated potentials, we find that there does NOT exist a 
DO£)*o (^d*0£)0^ molecular state if we use the experimen- 
tal value for the DD*?: coupling constant and a reason- 
able value around 1 GeV for the cutoff (A or ^/J3). The 
two approaches lead to the same conclusion. Bound state 
solutions with negative eigenvalues for the DD* system 
exist only with either a very large coupling constant (two 
times of experimental value) or a large cutoff (A 6 GeV 
or /3- 6 GcV^). 





P (GeV^) 


Eo (MeV) 


r-rms (fm) 


»"max (fm) 




= 0.59 


0.9 


-0.9 


2.0 


0.3 


y 


1.0 


-3.2 


1.2 


0.3 




= 0.7 


0.5 


-1.0 


2.0 


0.4 


a 


0.6 


-4.7 


1.0 


0.3 




= 0.8 


0.3 


-0.5 


2.9 


0.5 


9 


0.4 


-5.0 


1.0 


0.4 


9 


= 0.9 


0.2 


-0.4 


3.3 


0.6 


9 


= 1.0 


0.2 


-4.1 


1.1 


0.5 



TABLE VII: Solutions for various g and (5 in the case of smear- 
ing for the BB* system with OPEP. The lowest eigenvalues 
between -5.0 MeV and -0.1 MeV are selected. 





A (GeV) Eo (MeV) 


rrms (fm) 


r'niax (fm) 


g = 0.59 


2.5 


-0.5 


2.7 


0.3 


2.6 


-2.5 


1.2 


0.2 


5 = 0.7 


1.8 


-0.3 


3.8 


0.4 


1.9 


-1.9 


1.5 


0.3 




1.4 


-0.2 


4.4 


0.5 


3 = 0.8 


1.5 


-1.7 


1.6 


0.4 




1.6 


-4.9 


1.0 


0.3 


3 = 0.9 


1.2 


-1.1 


1.9 


0.5 


1.3 


-3.9 


1.1 


0.4 


3 = 1.0 


1.0 


-0.9 


2.1 


0.5 


1.1 


-3.7 


1.2 


0.4 



TABLE VIII: Solutions for various g and A in the case of 
FF for the BB* system with the total potential. The low- 
est eigenvalues between -5.0 MeV and -0.1 MeV are selected. 
Here po- = 0.76 is used. 



Because B mesons are much heavier, hence their ki- 
netic energy decreases which is helpful to the formation 
of the shallow BB bound state. In fact, our analysis indi- 
cates that there probably exists a loosely bound S-wave 
BB* molecular state. Once produced, such a molecu- 
lar state would be rather stable since its dominant decay 
mode is the radiative decay through B* — > B^. Experi- 
mental search of these states will be very interesting. 

In short summary, we have performed a dynamical cal- 
culation of the D^D*'^ system in the mature meson ex- 
change framework. Our analysis disfavors the interpre- 
tation of X(3872) as a loosely bound molecular state if 
we use the experimental coupling constant and a reason- 





P (GeV") 


Eo (MeV) 


rrms (fm) 


r-max (fm) 




= 0.59 


1.1 


-0.4 


2.9 


0.3 


9 


1.2 


-2.6 


1.2 


0.2 




= 0.7 


0.6 


-0.6 


2.6 


0.3 


9 


0.7 


-4.1 


1.0 


0.3 


9 


= 0.8 


0.4 


-1.4 


1.7 


0.4 


9 


= 0.9 


0.3 


-3.1 


1.2 


0.4 


9 


= 1.0 


0.2 


-1.8 


1.6 


0.5 



TABLE IX: Solutions for various g and /3 in the case of smear- 
ing for the BB* system with total potential. The lowest eigen- 
values between -5.0 MeV and -0.1 MeV are selected. Here 
Oct = 0.76 is used. 
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able cutoff around 1 GeV, which is the typical hadronic 
scale. Clearly more theoretical and experimental efforts 
are require to understand the underlying structure of the 
charming and mysterious X(3872) state. Maybe one need 
consider some more exotic schemes like the admixture of 
a cc charmonium and a DD* molecular state. Coupled 
channel effects will help further lower the energy of the 
system. 
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